Sequential Bifurcations in Sheared Annular Electroconvection 
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A séquence of bifurcations is studied in a one-dimensional pattern forming System subject to the 
variation of two expérimental control parameters: a dimensionless electrical forcing number 71 and a 
shear Reynolds number Re. The pattern is an azimuthally periodic array of traveling vortices with 
integer mode number m. Varying TZ and Re permits the passage through several codimension-two 
points. We find that the coefficients of the nonlinear terms in a generic Landau équation for the 
primary bifurcation are discontinuous at the codimension-two points. Further, we map the stability 
boundaries in the space of the two parameters by studying the subcritical secondary bifurcations in 
which m — > m + 1 when TZ is increased at constant Re. 



In typical nonlinear pattern- forming Systems, there 
is one easily adjusted quantity, the control parameter, 
whose variation results in a symmetry-breaking bifur- 
cation to an ordered pattern ||ï|. If a second control 
parameter is available, then it may be possible to tune 
the parameters to codimension-two points at which two 
modes with différent symmetries are equally unstable. 
Near such points, neighboring unstable modes compete 
with one another leading to low-lying secondary bifurca- 
tions even in the weakly nonlinear régime which is ac- 
cessible to perturbative theoretical treatment. There are 
few expérimental Systems whose patterns can be conve- 
niently manipulated by the independent variation of two 
control parameters. Examples include Rayleigh-Bénard 
convection (RBC) in a rotating cylindcr |^,||, convection 
in binary fluid mixtures 0,^, and Taylor Vortex Flow 
(TVF) with rotation of both the inner and outer cylin- 
ders In most cases, the eventual two-dimensionality 
of the patterns results in a very complex phenomenology. 

In this Letter, we study experimentally the bifurca- 
tions in a thin annular fluid film. The two-dimensional 
(2D) film is both sheared and driven into electroconvec- 
tion This System is unique in that it involves two 
easily adjusted control parameters, the applied voltage 
and the shear, with the bifurcations to uni-directionally 
traveling vortices only in one-dimension (ID). The annu- 
lar geometry makes the ID pattern effectively periodic 
and free of boundary efîects. Most other pattern-forming 
Systems are fluid-mechanically three-dimensional (3D) 
and the competing patterns often have différent symme- 
tries and dimensions. For example, in TVF there is a 
codimension-two point where Taylor and wavy vortices 
are simultaneously unstable. The former (latter) pattern 
is ID (2D). The greatly increased size of the parameter 
space made available by the third dimension makes it 
more difficult to examine experimentally the instabilities 
to and bctween the competing pattern states. Further, it 
is important to understand which physical processes dé- 
pend on the différent symmetries of the competing modes 
and which on their différent dimensions. In our System, 
symmetry alone détermines the nature of the bifurca- 
tions. 



The highly restricted geometry of our fluid film only 
allows the compétition of azimuthally periodic ID modes 
of varions m-fold symmetries. As a function of voltage 
and shear, we find a simple séquence of primary and 
secondary bifurcations. The primary bifurcation is the 
transition from the conducting state to the convecting 
State, whereas the secondary bifurcations connect one 
convecting state to another. This scénario can be mod- 
elled by coupled Landau amplitude équations which can 
be calculated completely by perturbative methods from 
the primitive electrohydrodynamic équations [0J^. The 
amplitude équations reduce to a single generic Landau 
équation when applied to the primary bifurcation. We 
experimentally demonstrate that the coefficients of the 
nonlinear terms in this équation vary discontinuously 
at codimension-two points. Our experiments also re- 
veal that the secondary bifurcations near codimension- 
two points are subcritical and suppressed by the shear. 

The expérimental System consiste of an annular, freely 
suspended film of smectic A liquid crystal. The layered 
smectic A phase is isotropic and newtonian for flows in 
the plane of the film, and maintains a uniform thickness 
even while flowing. Such a film is driven to electroconvect 
by a voltage applied between the inner and outer edges 
which sets up an unstable distribution of surface charge 
on the film's free surfaces |9j . A Couette shear is applied 
by mechanical rotation of the inner électrode. Previous 
studies of this instability explored the behavior of the 
primary bifurcation to electroconvection |^,^, confirm- 
ing that the weakly nonlinear régime of convection near 
onset is described by a complex amplitude équation of the 
Landau form for a single mode number m | pÔ[ . The con- 
centric électrodes have inner (outer) radii (ro), from 
which we define the radius ratio a = ri/ro. The film 
forms an annular sheet of width d = ro — Vi ^ 1mm and 
thickness s ~ O.l/im so that s/d ^ 10"'*. The fluid is 
treated as incompressible, strictly two-dimensional, and 
of uniform electrical conductivity. We dénote the fluid 
density, molecular viscosity, and conductivity by p, rj, 
and a. In a typical experiment, a steady circular Cou- 
ette shear was imposed by rotating the inner électrode at 
u! rad/s and a de voltage V was applied between the in- 
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ner and outer électrodes. Current-voltage characteristics 
were obtained for incrémental and décrémentai voltages. 
The state of the film is specified by the radius ratio and 
three additional dimensionless parameters. The control 
parameter TZ = e^V'^ / arjs^ , where éq is the permittivity 
of free space, is a measure of the external driving force. 
The Prandtl-like parameter V = e^rj/ pasd is the ratio of 
the time scales of electrical and viscous dissipation pro- 
cesses in the film. The Reynolds number Re = pcurid/r] 
indicates the strength of the applied shear and is the 
second control parameter. Détails of the expérimental 
apparatus and procédure can bc found in Rcfs. P,^-[rï[| . 
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FIG. 1. Plots of n versus e for a = 0.47. In (a) V = 15.3, 
Re = 0.94. In (b) V = 21.7, Re = 0.58. 

The primary bifurcation leads to a pattern of m vor- 
tex pairs. Onset occurs at a critical voltage V = Vc 
with the critical mode number m — rric becoming un- 
stable. Linear stability j|] predicts the values of rric 
and the critical control parameter TZc oc at which 
the film becomes marginally unstable. In gênerai, TZc 
and rric are functions of a, V, and Re. However, when 
Re = 0, TZc = TZc and nie ~ "^c are independent of 
T'. When Re > 0, it is found that Tic > ^c, for any 
a and T', Le. that the shear always acts to stabilize 
the conduction state and suppresses the onset of convec- 
tion. In addition for nonzero shear, nie < rric- Both 
thèse prédictions of linear theory have been experimen- 
tally confirmed. Our main expérimental tool is the to- 
tal electric current /, which is a measure of the rate of 
charge transport due to both conduction and convection. 
The reduced Nusselt number n is the ratio of the cur- 
rent due to electroconvection to that due to conduction 
and is defined a,s n — (I/Icond) — 1 = (I/cV) ~ 1- The 
conductance of the film c is obtained from the current- 
voltage characteristic of the film when it is not convect- 
ing. Similarly, the reduced control parameter e, is given 
by e = (TZ/TZc) - 1 = iylVcf - 1. The raw (y,/) data 
can thus be transformed to measurements of (e,n) [ p2[ . 
In thèse dimensionless terms, the primary bifurcation oc- 
curs at e = 0. When e is increased, the secondary bifurca- 



tions are transitions between flows with mode numbers 
TO ^ m -|- 1 Ip^ . Upon decreasing e, the transitions 
TO ^ m — 1 are observed. The mode changes can be 
observed directly, but rather crudely, using visual obser- 
vations of slightly nonuniform films, or films seeded with 
small particles. They appear very cleanly, if more in- 
directly, in the form of jumps in n as a function of e. 
Thèse observations are in qualitative agreement with an 
equivariant bifurcation theory Increasing (decreas- 
ing) e results in secondary bifurcations at ei < 62 < . . . 
(e_i < e_2 < . . .) corresponding to mode transitions 
TO^TO-|-1— (. ..-^to4-1— > to). 

In Figure ^ are illustrated two représentative examples 
of secondary bifurcations. Figure |](a) shows (e,n) data 
for a sheared film with subcritical primary and secondary 
bifurcations. The secondary bifurcation, which results in 
one additional traveling vortex pair, appears at ei 0.22 
for increasing e. When e is decreased, the removal of 
the vortex pair occurs at e_i 0.12. Convection al- 
together ceases at e ~ —0.1. Figure |^(b) shows a case 
where, for increasing e, the primary and secondary bifur- 
cations appear as a single strongly subcritical bifurcation 
at e ~ 0. However, for decreasing e, e„i is distinct from 
e ^ —0.11 where convection stops. That shear sélects 
the mode which is linearly unstable is well established 
in the case of the primary bifurcation It has been 
shown from linear stability analysis that the marginally 
unstable mode decreases from to|Î — 1 ^ to^ — 2 . . . 

as Re is increased. Hence, for fixed a and T', there are 
codimension-two points at particular Reynolds numbers, 
Rem.m-i, at which both to and m — 1 are marginally 
unstable |ï^. For Re < Kem,m-i, the primary bifurca- 
tion is to a convecting state with mode number m. The 
first secondary bifurcation occurs at e = ei > to mode 
number to -|- 1. For Re > Rem,„i-i, the primary bifur- 
cation is to mode m — 1 . The first secondary bifurcation 
is at e = Êi > to mode to. Clearly, it foUows that as 
Re — ^ Rem,m-i from above, ei — > 0, resulting in the sit- 
uation depicted in Fig. |ï|(b). Our expérimental protocol 
was always to hold Re constant and vary e. Neither ex- 
periment nor theory has addressed the question of what 
mode transitions occur if Re is varied for constant e. 

Figures ^(a)-(h) show a séquence of (e, n) data plots 
for a = 0.56, T' = 76, and a range of closely spaced 
values of Re. In Fig. ||(a), the data at Re = 0.124 
shows a subcritical primary bifurcation with mode num- 
ber TO = 8. For the range of e investigated, secondary 
bifurcations were not encountered. In Fig. ^(b) , the data 
at Re = 0.142 shows a subcritical primary bifurcation 
with m = 7 foUowed by a subcritical secondary bifurca- 
tion (m = 7 8) at ei 0.26. Hence at some Reynolds 
number, 0.124 < Regj < 0.142, there is nascent at 
e = the mode change 8^7. At Re = 0.160, the sec- 
ondary bifurcation is suppressed to higher ei ~ 0.48, as 
shown in Fig. ^c). The séquence of events in Figs. ||(a) 
through (c) repeats with increasing Re. InFig. |(f), data 
at Re = 0.214 shows a primary bifurcation to to = 6 and 
a secondary bifurcation at ei ~ 0.07 to m = 7. Hence a 
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second codimension-two point at 0.196 < Rey.g < 0.214 is 
traversed. At Re = 0.231 and Re = 0.249, the secondary 
bifurcation is suppressed to ei ~ 0.11 and ei ~ 0.16 re- 
spectively, as is shown in Figs. ||(g) and (h). The data 
obtained while e is being reduced has equally interest- 
ing behavior. For example, for increasing e in Fig. |2|(f) 
there is a m = 6 — * 7 bifurcation, but for decreasing e 
the m = 7 convecting state exits directly to conduction 
by a single discontinuons transition, without recourse to 
the TO = 6 state. However, for the slightly higher Re in 
Figs. ^(g) and (h), the to = 6 state is revisited while e 
decreases. 
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FIG. 2. A Re-sequence demonstrating mode changing bi- 
furcations. Plotted are n versus e at a = 0.56 and V = 76 for 
a séquence of increasing Re. The open (filled) triangles are 
for increasing (decreasing) e data. 

For nonzero shear, symmetry and perturbative analy- 
ses of the primitive équations lead to the amplitude équa- 
tion 

T(Am ^ iO'ImA^ri) = e(l + iCOm)Am 
-50m(l + iC2m)\Am\ A„i - C/lm(l + )\An\'A 
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and its counterpart under m ^ n symmetry In the 
above, ail of the coefficients are real. The steady state 
solutions of the real part of Eq. 3.1 are consistent with the 
expérimental observations of a state m being replaced by 
a state n = to ± 1 via a discontinuons transition Q . To 
further interpret the data in Fig. ||and better understand 
the behavior of the primary bifurcation as the Reynolds 
number is increased through Rcsj and Re7_6, we consider 
the real and steady state part 
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of Eq. |0.l| , where A is real, g = gom, and h = hom- We 
have truncated at quintic order in the expansion. Since 
the reduced Nusselt number n = we can dé- 

termine g and h for each primary bifurcation in the sé- 
quence in Fig. H. Our results are plotted in Fig. |. We 
find that both g and h are strongly discontinuons as Re is 
increased through the codimension-two points Resj and 
Rc7^6- The discontinuity in both g and h is smaller at 
Reg_7 than at Rey.g. This raises the possibility that thèse 
discontinuities may vanish in the limit m — > cx) and so are 
a feature of the finite size and discrète mode structure of 
this System. 




FIG. 3. The coefficients g and h in Eq. D.2 versus Re for the 
data in Fig. 2. The dashed lines locate the codimension-two 
points. Solid lines are linear fits to the data. 

A minimal description of a state consists of identifying 
the value of its mode number m. AUowing for hystérésis, 
each state persists over certain overlapping ranges of a, 
■p, Tt, and Re. Here, we focus on fixed a and V for 
In Fig. 



simplicity |15|. In Fig. U we map the solutions in the 
subspace defined by (e,Re, a — 0.56, 'P = 76), where the 
reduced variable e is used for convenience. Figure |^ shows 
the values of e at which conduction and convection with 
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various m were found, for given values of Re. Below 
(above) the lowermost dashed Unes (horizontal solid line), 
only conduction (convection) occurs. Between thèse lines 
both conduction and convection are observed, depending 
hysteretically on whethcr e is increasing or decreasing. A 
tricritical point(TC) is shown at Re ^ 0.11. For Re 
0.11 the onset of convection is supercritical and hence 
non-hysteretic [0. 
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ondary bifurcations appear near codimension-two points. 
Using a generic Landau model, we have quantitatively 
analyzed the behavior of the primary bifurcations on ei- 
ther side of thèse points. We have found that the co- 
efficients of the cubic and quintic terms of the Landau 
équation are strongly discontinuons at the codimension- 
two points. We expect that this scénario occurs over large 
régions of our parameter space and may well be a gên- 
erai property of discrete-mode pattern-forming Systems 
such as TVF. ^ Unlike TVF however, this System only 
has bifurcations between one-dimensional patterns with 
simple m-fold symmetries. The parameter space is char- 
acterized by branches that demarcate différent stable and 
multistable domains. The surprizingly intricate bifurca- 
tion structure of this highly geometrically constrained 
instability is amenable to detailed quantitative analysis 
starting from the underlying electrohydrodynamic équa- 
tions, [pj As such, it is an idéal venue for expérimental 
and theoretical studies of bifurcations in finite nonlinear 
pattern-forming Systems. 

We thank R. E. Ecke for helpful discussions and con- 
structive comments. This research was supported by the 
Canadian NSERC and the U.S. DOE (W-7405-ENG-36). 



FIG. 4. Plot of the parameter space of patterns with mode 
numbers m for a = 0.56 and V = 76. Open (closed) symbols 
indicate conduction m{m — > conduction) and m — > m -|- 1 
(m -f 1 — > m) transitions encountered for increasing (decreas- 
ing) e. The numbers are the observed m. The région below 
(above) e = is multistable with conduction (m and m + 1 
States). TC is a tricritical point that connects the super- 
critical and subcritical branches. CoD2 are codimension-two 
points. The solid lines for e > are the secondary stability 
boundaries across which m — > m + 1. Each of the dashed lines 
for e < connects the positions of the saddle-nodes for mode 
m. 

Within the convecting régime, the boundaries of the 
secondary bifurcations between m = 8, 7 and m = 7, 6 
show upwards trends indicating that an analogous sup- 
pression occurs for the modes m = 8 and m = 7, even 
after they cease to be the primary mode of bifurcation 
. The suppression of the secondary modes increases 
more rapidly with Re than when they were the primary 
mode, with larger m modes suppressed more rapidly than 
smaller m modes. Inside the hystérésis loop of the pri- 
mary bifurcation, we found as many as three possible 
distinct states, one being conduction. (See Figs. ||(b), 
(c), (g), and (h)). The entire portrait presented in Fig. ^ 
dépends perhaps strongly on a and weakly on V and is 
likely to be more branched (contain more transitions) for 
larger a. It should be possible to generate theoretically 
such a tree from the primitive é quat ions, once they have 



been reduced to the form of Eq. 3.1 



We have shown that in our discrete-mode System, sec- 
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